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Abstract 

We introduce and study a new discrete basis of gravity constraints by making use 
of harmonic expansion for closed cosmological models. The full set of constraints 
is splitted into area-preserving spatial diffeomorphisms, forming closed subalgebra, 
and Virasoro-like generators. Operatorial Hamiltonian BFV-BRST quantization is 
performed in the framework of perturbative expansion in the dimensionless param- 
eter which is a positive power of the ratio of Planckian volume to the volume of the 
Universe. For the {N +\) - dimensional generalization of stationary closed Bianchi-/ 
cosmology the nilpotency condition for the BRST operator is examined in the first 
quantum approximation. It turns out, that certain relationship between dimension- 
ality of the space and the spectrum of matter fields emerges from the requirement 
of quantum consistency of the model. 

PACS numbers: 04.60.Ds; 98.80.H; 04.60.-m 

Keywords: quantum cosmology, constrained dynamics, BFV-BRST method. 



1 Introduction 



At the present time we observe the explosive development in quantum gravity and string 
theory which have many-sided interconnections [|I|. From the standpoint of the con- 
strained dynamics p[ these theories manifest a certain similarity originating from the 
reparametrization invariance which is a gauge symmetry of the both. Meanwhile, the 
theories are usually treated in an essentially different way in respect to the perturbative 
expansions. Various well-advanced string models, having a trend to be thought about as 
dual ones, are considered nowadays as expansions around different limiting values of a 
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coupling constant of an unknown fundamental superstring theory 0] . Oppositely, gravity 
has a widely recognized candidate for a fundamental theory that is Einstein action or, 
perhaps, its some or another extension, but the perturbative treatment is a longstanding 
unsolved problem of the General Relativity. 

The main intention of the present work is to undertake an attempt of the Hamilto- 
nian BFV-BRST quantization of closed cosmological models by exploiting the keenship 
between the string and gravity gauge symmetries. 

Apparently, a phase-space constrained description is very similar for string theory and 
gravity while they are being carried out locally with dynamical variables and constraints 
defined at a point of a space-like section of world-sheet or spacetime, respectively. How- 
ever, there is an essential distinction in habitual approaches to quantization of strings and 
gravity. 

Quantizing strings, one may work with a discrete set of modes, describing string 
excitations, and a discrete set of constraints, which form Virasoro algebra. In doing so, 
one observes that quantum commutators of the constraints acquire central extension, and 
the quantum consistency of the theory is provided for the critical values of dimensionality 
and intercept. One of the most efficient methods to observe this phenomenon is the BFV- 
BRST Hamiltonian quantization of constrained systems ^, ^, 0. In the framework 
of this method one can show how the nilpotency condition for the BRST charge defines 
the values of critical parameters for strings p, 0, membranes |]10[ and W3 gravity |11 



Actually, the use of discrete set of modes and constraints is merely a convenient tool for 
taking into account the topology of the compact world-sheet spacelike section, when the 
method is applied. 

Canonical quantization of gravity and cosmology is usually treated by either perturba- 
tive consideration of path integral or operatorial quantization according to Dirac scheme 
with constraints imposed onto physical states. Basically, such a treatment implements 
continuous basis of dynamical modes and constraints, which can satisfactory describe the 
dynamics locally but are hardly well-defined globally on compact spacelike sections of 
the spacetime. Meanwhile, the consideration of cosmological perturbations in terms of 
globally defined harmonics was used for the investigation of different problems in classical 
and quantum |T3|, cosmology. However, a globally defined discrete constraint basis 
for quantum cosmology of closed Universes has not been hitherto discussed. Consistent 
canonical quantization requires both the dynamical modes and constraints to be globally 
defined, since this is essential for the adjustment of the ordering of the modes to the 
quantum algebra of the constraints. For example, in the case of the Wick ordering of 
bosonic string modes, the quantum corrections to the constraint commutators give rise 
to the central extension of the Virasoro algebra representation in the Fock space. Then 
the quantum constraint algebra can be made consistent provided the proper basis of con- 
straints is chosen (L and L are employed instead of habitual for cosmologists lapse and 
shift generators H±, H\\ and moreover, only one half of L„ and Lm should annihilate 
physical states). At the mean time, the straightforward pursuing the Dirac scheme, with 
i^x and H\\ being imposed on the physical subspace, leads to a contradiction. 

In this paper, we suggest a new scheme for canonical quantization of the closed cos- 
mological models. The distinguishing features of this treatment are as follows: the de- 
composition of both dynamical variables and constraints into harmonics which are the 
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eigenfunctions of the Laplace operator for the maximally symmetric space of given topol- 
ogy; separation of subalgebra of area-preserving diffeomorphisms from a total set of gravity 
constraints whereas the rest of the constraints form Virasoro-like generators; quantiza- 
tion of this constrained theory is performed in the framework of the Hamiltonian BFV- 
BRST formalism with due regard to the ordering of gravity, matter and ghost harmonics; 
studying the quantum nilpotency condition for the BRST operator, we apply a certain 
perturbative expansion of the constraints and structure functions with a small parameter 
Ip/V^^^, where Ip is a Planck length and is a spatial volume of the Universe, is a 
dimensionality of the space. Implementing this scheme, we focus the consideration mostly 
at the stationary Bianchi-J type Universe, although any other closed cosmological model 
could be treated in a similar way. The specificity of a particular topology is encoded in the 
constraint algebra structure functions which appear to be expressed via Clebsch-Gordan 
coefficients of the corresponding symmetry group representation. In the case of Bianchi-J 
cosmology, these coefficients are the simplest, as the group is f/(l)^. Although we take 
this particular type of a closed cosmological model for the sake of technical simplicity. 



it is, however, of a certain physical interest |15|. On the other hand, A^-torus allows to 



exploit a straightforward analogy with the closed bosonic string sigma-model when the 
nilpotency of the quantum BRST charge is examined. The most striking outcome of this 
quantization procedure, being applied to the Einstein gravity coupled to the matter fields, 
is that the number of the matter degrees of freedom appears to be correlated with the 
dimensionality of the space A^. For example, when there are d mass less scalar fields only, 
already the first quantum correction gives rise to the relation between d and A^: 

d = 30 + ^{N + l){N-2), (1.1) 

which is a necessary condition for nilpotency of the quantum BRST charge. Mention the 
curious fact: if one puts A^ = 1 (thereby torus reduces to a circle and the constraints form 
the Virasoro algebra) then relation (1.1) results in d = 25. This result can be naturally 
understood from the standpoint of the string sigma-model without Weyl invariance |1^ . 



As a matter of fact, bosonic string theory should not be thought about as one-dimensional 
limit of the Einstein gravity coupled to a set of d scalar fields, because string possesses 
an extra gauge symmetry, namely Weyl invariance, besides the diffeomorphisms, whereas 
the mentioned a-model does not have Weyl invariance. As is known |jl6[ the critical 



dimensionality of the a-model is 25, in contrast with 26 for strings. The discrepancy 
in these dimensionalities can be easily explained if one remembers that Weyl invariance 
gauges out conformal mode of 1+1-metric, while in the non-conformal cr-model case this 
mode contributes to the Virasoro generators on an equal footing with string excitations 
16| . Thus, in some sense, the critical dimensionality for the both cases is 26, but for the 



(T-model this number consists of 25 string coordinates and 1 gravity conformal mode. 

In its basic features the suggested scheme should, probably, be actual for other closed 
cosmological models, including those with non- stationary classical background. But the 
critical relations analogous to (1.1) may change their form for each particular model. 

The structure of the paper is as follows: in the second section we describe a general 
decomposition of the gravity constraints into the discrete set of area-preserving diffeomor- 
phisms and Virasoro-like generators and sketch the way of the representing the structure 
constants of the constraint algebra via Clebsch-Gordan coefficients of the corresponding 
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symmetry group; in the third section we suggest the spht of integer vectors on A^-torus 
into two classes which furnishes a means of choosing the Wick ordering for a proper part 
of dynamical modes and constraints; in sec. 4 we write down explicitly the decomposition 
of gravity constraints in new discrete basis and the corresponding structure constants for 
closed A^-dimensional Bianchi-J cosmological model (A^-torus); in sec. 5 we define per- 
turbative expansion of the constraint operators and calculate first quantum corrections to 
the constrained algebra; sec. 6 is devoted to the investigation of the structure of BRST 
operator and examination of its nilpotency condition, as a result, we come to the critical 
relation between the dimensionality of the space and the spectrum of matter fields; sec. 
7 contains conclusions. 



2 Gravity constraint algebra for a closed Universe. 
Harmonic expansion, 
area-preserving difFeomorphisms 
and Virasoro-like generators 

Hamiltonian description of any reparametrization-invariant theory with metric, being 
considered as a dynamical variable, contains a set of first-class phase-space constraints 
if_L and Hi, i = 1, . . . , N where is a dimensionality of space. H±, usually called super- 
Hamiltonian or Wheeler-DeWitt operator (the canonical Hamiltonian vanishes on the 
constraint surface for this case), generates lapse transformations of spacelike hypersurface, 
while Hi, called super momentum, generates spatial diffeomorphisms [|17|. If the gravity 
action does not contain higher derivatives of metric, these constraints can be represented 
in the following form 

H± = /p~"'^G'jj_jt;7r*%^' — j^_^ gR + H± matter, (2.1) 

Ip 

Hi 2,gijTC'^ {^Qik,m ~\~ 9im,k Qkm,i)'^ ~\~ Hi matter i l^'^) 

where gij is a metric on a spacelike section of spacetime, vr*-^ is a momentum conjugated to 
metric and representing extrinsic curvature of the spacelike hypersurface, Ip is the Planck 
length. The DeWitt supermetric Gij^ki has the following form: 

Gij^ki = - (^gikQji + giiQjk — _i dijOki^ , (2.3) 

while its inverse looks as 

H± matter and Hi matter are the matter field contributions to super-Hamiltonian and 
supermomenta, respectively. Their explicit form depends upon a specific set of the matter 
fields and their interactions, although the involution relations between the constraints are 
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not effected by the structure of matter part of tliese expressions, because tfiese relations 



liave a purely geometric origin [IS]. Tliese involution relations are 



{H^{x),H^{x')} = gix)g'\x)H,ix)^5{x,x') 
-g{x')g'^{x')H,{x')^6{x,x'), (2.5) 

{H,{x),H^{x')} = H^{x)-^^S{x,x') 

-H^{x')^^5{x,x'), (2.6) 
{H,{x),H,{x')} = H,{x)^6{x,x') 

-H,{x')^6{x,x'). (2.7) 

Here {, } denotes conventional Poisson brackets. The expression for super-Hamiltonian 
( p.l|) differs from the conventional one by the factor g^^'^, thereby the involution relations 
( p.5|) , (|2.6|) also have distinctions from the commonly used. Such an unusual definition of 
super-Hamiltonian (which, however, sometimes occurs |l^) is convenient for our purposes. 
Notice, that the involution relations (p.5|)- ( |2.7| ) could not be reproduced without gravity 
dynamical variables, because in its essence they express the gravity gauge symmetry 
algebra in phase space [pO| ]. 

Considering compact manifolds, one can express phase space variables and constraints 
among them via a discrete set of coefficients of harmonic expansion. The functional ba- 
sis of the expansion is formed by eigenfunctions of the Laplace operators defined in the 
maximally symmetric space of a topology under consideration. For the case of BianchiJ 
cosmology these functions are simply elements of Fourier expansion on a torus, for the 
case of Bianchi/X cosmology - hyperspherical functions ||12[ . From the viewpoint of lo- 



cal consideration this harmonic representation is completely equivalent to the original 
continuous parameterization of the phase space and constraints, but the benefit of this 
treatment is that the expansion, being well-defined on the manifold, reflects explicitly 
all the topological distinctions between different types of cosmological models. The sig- 
nificance of this topological specificity becomes obvious at comparing open and closed 
bosonic strings. Locally, both the models are described by the same canonical variables 
and constraints, whereas the difference transparently manifests itself in the spectrum of 
modes and constraints upon the Fourier expansion. Really, the dynamical modes of closed 
string are left- and right-moving excitations aj^, am subjected to two copies of Virasoro 
constraints L„ and L^, while in the case of the open string the dynamical modes are the 
standing waves with a single set of Virasoro constraints L„. 

Now, we shall try to push the consideration of the constrained Hamiltonian formalism 
along the lines of harmonic representation for the closed cosmological models. The super- 
Hamiltonian can be written as 

H^{x)=J2hPq^''\x), (2.8) 

(n) 
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where Q^"'\x) is a scalar-type harmonic of the corresponding Laplace operator, labeled 
by a multi- index (n) . The discrete set of constraints H^^^ is given by 

= j dVH^{x)Q<^\x), (2.9) 

where dV denotes the covariant volume element of the manifold; the harmonics are nor- 
malized 

j Q(")(x)Q*(")(x) = (2.10) 
The structure of supermomentum expansion is more complicated: 

Hi{x) ^Y.H\\^da''\x) + Y.H^^^S\^\^) > (2-11) 

(") {A) 

here 5'^-'^'' [x) stand for the transverse vector harmonics obeying the relations 

y^S\^\x)=Q. (2.12) 

So, two types supermomentum constraint harmonics appear: the longitudinal H\\ and the 
transversal if^^^ labeled by another set of multi-indices {A). The inverse relations for H\\ 
and are 

i/{") = j dV H,{x)^^^^^ , (2.13) 

ii-(^) = j dV Hi(x)S*'^^\x) , (2.14) 

where A*^"^ is an eigenvalue of the Laplace operator corresponding to the eigenfunction 
(^"•"^(a;). Transversal vector harmonics are normahzed 

jdVS\^^S'^^^^8^^^^^\ (2.15) 

Transversal constraint's harmonics if^^^ define corresponding continuous subset of super- 
momentum constraints 

^,(x) = E//(^)5f^)(x). (2.16) 

(A) 

Constraints Hi{x) are obviously divergenceless 

V^^i(x) = 0, (2.17) 

where the covariant derivative is defined with respect to the metric of the maximally 
symmetric space, which we shall call as a background metrics in what follows. 

It is easy to check that the Poisson bracket of Hi{x) among themselves gives Hi{x) 
again 

{H,{x) , H^{x')} = H^{x)^S{x,x') - H,{x')^5{x,x') (2.18) 

Usually, it is suggestive to consider the contractions of constraints to arbitrary "weight- 
functions". In the case of Hi{x) , these constraint functional 

j dVH,{x)r{x)=H{f), (2.19) 



where can be chosen divergenceless, form the closed algebra with respect to Poisson 
bracket 

mh) , HU2)} = Hiifi, f2]) , [/i, = fi^n - fi^fi ■ (2.20) 

The Lie bracket of divergenceless vector fields /i,/2, apparently, is a divergenceless too. 
The commutation relations (|2.20|) are typical for generators of area-preserving diffeomor- 



phisms introduced in the paper ||2T[ and being intensively studied in relationship with 



p-brane theories [^. However, the transformations generated by H do not preserve the 
volume of the physical space, i.e. {H,g} 7^ 0, where g = detgij is a genuine (not a 
background) metric. In principle, one can introduce divergenceless projection of super- 
momentum with respect to the genuine metric gij, and such a projection does generate 
area-preserving transformations. However, these projections of the constraints do not 
form a closed algebra with respect to Poisson bracket. This disclosure is caused by the 
fact that the projector depends upon the metric and so it is not invariant with respect to a 
metric shift. So, the Poisson bracket of two area-preserving supermomentum constraints 
does not give an area-preserving projection of a supermomentum constraint again. For 
our purposes it is essential to employ the constraint basis, with the closed subalgebra of 
spatial transversal diffeomorphisms, and we adopt for this reason the basis of longitudinal 
and transversal constraints. The latter we shall refer to as area-preserving, because they 
form the closed subalgebra of this type, although they employ a background metric in- 
stead of a genuine one. The longitudinal part of supermomentum will be combined with 
the super-Hamiltonian discrete constraints (|2.9| ) to construct Virasoro-like generators. 
In continuous basis these generators read 

H^{x) ± ^V^/7,(x), (2.21) 

where the covariant derivative V* and Laplace operator A are defined in respect to back- 
ground metric. Similarly, area-preserving diffeomorphisms in continuous parameterization 
are 

1 



5}-V,-V^ji/,. (2.22) 

Dynamical variables are also expanded into the appropriate set of harmonics. For 
example, a scalar field expansion is 

^(t,x) = $:^(")(t)Q(")(a;), (2.23) 

(n) 

where Q*^"-* is the same set of scalar harmonics as in the expansion of H±{x) ( |2.8| ). The 
expansion of metrics is more complicated: 

g,,{t,x) =Y.a^-\t)g§\x)Q^-\x) + T.b^-\t) - ^ Q^-\x) 

(n) (n) \ / 

+ Y.<^^^\t){v4^\x) + V,S\^\x)) + j:S^\t)Gi;\x), (2.24) 

(A) (a) 

where glj\x) is the background metric, G\j\x) is the set of transverse traceless tensor 
harmonics: 

. = 0; V'Gij = 0. (2.25) 
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Similar expansions can be written for the rest of canonical variables. For the case of 



this multipole expansion (|2.24|) was described in detail in |12 



The expansion for Qij (|2.24| ) is involved not only into super-Hamiltonian (|2.1|) and 



supermomenta ( |2.2| ), but also into the structure functions of constraint algebra via relation 
( p.5|) , thereby the multipole modes a^'^\t) , b^"'\t) , c^^\t) and S"\t) enter the involution 
coefficients. However, in the framework of the perturbative consideration adopted in this 
paper for the calculation of quantum corrections, the first approximation can be obtained 
neglecting the multipoles from the involution structure functions. This will be explained 
in more detail in Sec. 5, where quantum BRST charge will be studied for the case of 
Bianchi-J cosmology. 

Let us discuss the multipole independent part of the structure functions. One may 
easily see that they are proportional to the integrals of three harmonic functions of such 
types: 

J dVQ^''\x)Q^"'\x)Q*^^\x), 
dVQ^''\x)Q^"'\x)Sf^\x), 
dVQ^''\x)sf^\x)sf''\x), 

dVSi^\x)S^f\x)S:^''\x). (2.26) 



The integrals (^^261) can be expressed via bilinear combinations of Clebsch-Gordan 
coefficients of the corresponding symmetry group by making use of Wigner-Eckart theorem 
(see e.g. H). 

If the expansion of constraints into the harmonics of dynamical modes is considered, 
then the coefficients at bilinear combinations of multipoles are expressed in terms of the 
similar integrals, which in their turn can be reduced to the Clebsch-Gordan coefficients. 
As a matter of fact, only quadratic combinations of harmonics in the constraints may 
contribute to the constant part of the first quantum correction to the constraint involution 
relations. This contribution can usually be very instructive, e.g.: for the case of string 
P, ^, membranes [|lOl or W^-grayitj [|lT| it is the correction, which fixes the critical 



parameters. 

From the viewpoint of obtaining the critical parameters for closed cosmological mod- 
els, the method is basically the same for any group, but the simplest case of Bianchi-J 
cosmology allows to avoid the cumbersome expressions involving summing of bilinear 
combinations of Clebsch-Gordan coefficients. 



3 Wick ordering and basis of integer vectors on A^- 
torus 

In contrast to ordinary quantum mechanics, quantum theories with an infinite number 
of degrees of freedom may be inequivalent for a different choices of operator symbols. Of 
course, different symbols of operators (different orderings) are formally linked by unitary 
transformations, but an appearance of ultraviolet divergences may destroy the unitarity 
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of the transformation. For example, Wick symbols of the Virasoro constraints in bosonic 
string theory yield the nonvanishing finite central extension, whereas Weyl symbol of the 
same constraints does not reveal this correction to the commutation relations. So, the 
choice of ordering appears to be not technical, but fundamental step, defining the theory. 

The common wisdom is that the modes of oscillating behaviour should be quantized 
by Wick ordering, while those, having non-oscillating dynamics (usually, zero modes) are 
ordered by Weyl rule. 

To provide an opportunity of Wick quantization of the part of the constraints and 
Fock representation for the model we should introduce some kind of split of the set of 
A^-dimensional wave vectors on A^-torus. In the case of string the similar problem can 
be resolved in a very simple way because there are one-dimensional vectors (i.e. integer 
numbers) n. One merely considers positive integer numbers as the numbers corresponding 
to the annihilation operators a„ and creation operators and negative integer numbers 
as numbers corresponding to the creation operators a+ and annihilation operators an- 
In turn, ghost operators corresponding to Virasoro constraints L„ and where n 

is positive are treated as creation operators while the ghost operators corresponding to 
constraints and L„ are treated as annihilation operators. 

An attempt to consider Virasoro-like constraint subset for multidimensional cases 
should imply some splitting of integer vectors 

n = (ni,n2, ■ ■ ■ ,nAr) (3.1) 

into positive and negative ones. The plausible way to do it, from our point of view, is 
suggested in the paper |^ . 



Let us fix in some way the order of main axis of our torus. Then, if the first wave 
number rii is positive, we shall call the vector positive and if rii is negative we shall call 
the vector negative. If rii = we should turn to the second integer number n2 and to 
call the vector n positive or negative depending on sign of n2. If both numbers ni and 
n2 are equal to zero, the sign of vector n is determined by the sign of and so on. One 
can easily see that all the vectors n must belong to one of these two classes excepting the 
zero- vector (0, 0, ■ ■ ■ , 0). Below we shall implement a notation 

sgn{n) = 1 , if n is positive 

sgn(n) = —1 , if n is negative. (3.2) 
As a matter of fact we shall consider the integer vectors weighted by radii of A^-torus: 

n=r-^,...,T)- (3.3) 



Ri R 



■N 



In what follows we shall basically use the so called "last element" limit, i.e. we shall 
reduce our calculations mainly to the case then only number un differs from zero. It will 
give us an opportunity to make all the calculations explicitly. 
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4 Constraint algebra of Bianchi-/ model 
in harmonic representation 

In this section, we write down explicitly Virasoro-like and area-preserving diffeomorphism 
constraints for Bianchi-7 model in harmonic representation and study their commutation 
relations. 

i^(g) = ^(gx(g) + ^^"'">;^-'"' ), (4.1) 
m = i U(n) - "^"'"1^1^''"' ) , (4.2) 

where 

H^{n) = j d^xH^{x)exp{inx), (4.3) 

Hi(n) = J d^xHi(x)ex.p(inx), (4.4) 

Here, positive and negative constraint harmonics L{n),L{n) and L{—n),L{—n) are 
complex conjugated: 

L*{n) = L{-n), L*{n) = L{-n). (4.5) 

Let us write down the involution relations among these constraints, omitting the grav- 
ity harmonic contributions from structure coefficients: 

im, L(m)] = UtS^Lin + m) + ui$;,fi^L{n + m) + C/^^SSJ^.(n + m) 

1 / ^ sgn{n)sgn{m)rnn\ 
4 y |m||n| J 

X < lsgn[n) \n\ — sgn[m)\m\ H ^ ^ ^ ^ ^ I L[n + m) 

( SQTi(tI ~\~ TTl] (tI^ — ?72^) \ — 

-|- I sgn{n)\n\ — sgn{rh)\rh\ p — L{n + rh) 







\n H 


- m 



+ [ni-mi ^ ]Hi{n + m)}, (4.6) 

\ [n + my J J 



[Lin), Lim)] = Ut^^Lin + rn) + t/fgt'^i)L(n + m) + U^i^^H.in + rh) 

1 / sgn{n)sgn{rh)rhn\ 
4 y |m||n| J 

if /^M^i , /-M-i , sgn{n + m){ff -m^)\ ^ 
X < I sgn[n)\n\ + sgn[m)\m\ H — — L[n + m) 







\n 4- 


- m 



-I- sgn{n)\n\ + sgn{m) \m\ — — L{n -\- m) 







\n H 





rii - m, - ^ ^ Hi{n + mU , 4.7 

(n + m)^ J J 
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[L(n), Urn)] = Ui^;yll^L(n + rn) + ull^L{n + rn) + U^^^^ + -) 
1 / sgn{n)sgn{m)mn\ 
4 y |m||n| j 

if sgn(n + m)(ff - m^)\ ^ ^, 

X < sgn[mj\m\ — sgn[n)\n\ H — — L[n + m) 

[\ \n + m\ J 

( /^M^i /^M^i sgnin + m)iff - m?')\ - . ^ ^, 

+ sgn[m)\m\ + sgn[n)\n\ — — Lin + m) 

\ \n + m] J 

Cin^ ~ fn?)ini + mi)\ ^,1 , 

Ui-nii-^- \' 'l\Hi{n + m)\ . 4.8 
[n + my j J 

In the right-hand side of the relations ( [4.6|) - ( [4.8|) appear supermomenta constraints 
ifj(n + m) contracted with vector 

(n^— m^)(n + m) 

n-m ^ , 

[n + my 

which is orthogonal to (n + m). Thus, we have seen, that in the involution relations for 
Virasoro-like constraints L(n) and Lin) appear area-preserving diffeomorphisms. So, the 
constraints L, L do not form a closed algebra. 

It is convenient to describe area-preserving diffeomorphisms as 

E^(n) = ViHi{n), vn = 0. (4.9) 

introducing an arbitrary vector v, being orthogonal to the integer vector n. The other 
involution relations are as follows: 



m] 



>■ "\ V 'i H(v,n)L(m) ^ ' H(v,n)L(m) ^ ' H(v,n)L{m) ^ 

I f sgn{m)sgn{n + m)\m\\ ^ 

= — (mv) 1 H — L{n + m) 

2^ ^\ \m + n\ J ^ ' 

1 f sgn{m)sgn{n + m)\'m\\ - ^ 
— [mv) 1 — — L[n + m) 

2 \ \m + n\ J 

Isqn(m) f,^^. rn^(mv)(ni + mi)\ ,^ ^, 

+ - I ^1 ' {rhn)v, - {mv)mi + \^'y^., ' Hi{n + m), (4.10 

2 \m\ \ [n + m)"^ I 



[H,{n), L(m)] = U'^ftfl .Mn + m) + ulft"^). Mn + rh) + U^fjy.f]Mn + m) 

L \ /J H(v,n)L(m) ^ ' H(v,n)L(m) ^ ' H{v,n)L{m) ^ ' 

1 f sgn{m)sgn{n + m)\m\\ ^ 
= — [mv) 1 — — L[n + m) 

2 \ \m + n\ J 

I f sgn[m)sgn[n + m)\m\\ - ^ 

— (mv) 1 — : — L{n + m) 

2^ '\ \m + n\ J ^ ' 

1 sqn[m) f,^^. /^^x rh"^ [mv) [rii + mi)\ ^^.^ ^, 

+ - I ^1 ^ {mn)v, - [mv)mi \^ Hi[n + m), (4.11) 

2 \m\ \ (n + mj^ J 
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{mw)H^{n + ffi) 

Concluding this section we should notice that the formulae ( |4.1| ) and ( |4.2| ) contain an 
ambiguity concerning zero modes i^(0) and 1^(0) which have to be described additionally. 
It is convenient define them as 

m = \{H^{Q)+H^m (4.13) 

and 

L(0) = l(i/^(0)-i7^.(0)) (4.14) 

to provide a proper correspondence with a genuine Virasoro algebra in the limiting case 
N =1. 



H{u,n+rn) 
H(v,vecn)H 

nv)HJn + m). (4.12) 



H(v,vecn)H(w,vecm) ^ ' 



5 Quantum corrections to the algebra 
of Virasoro-like constraints 

In Sec. 3 we have introduced the split of integer vectors on A^-torus into positive and 
negative ones. Now we shall apply this decomposition to define creation and annihilation 
operators for gravity and matter fields, and to split the Virasoro-like constraints into 
positive and negative ones. The latter decomposition will further define the natural Wick 
ordering for the ghost variables upon the construction of quantum BRST charge. Then we 
calculate the c-number part of the quantum correction to the commutator [L{n), L{—n)], 
which is of crucial importance from the standpoint of the nilpotency condition on BRST 
operator. 

To begin with, we consider the massless scalar field contributions to constraints Hj_{n) 
and Hi{n) which, in turn, define Virasoro-like constraints. Scalar contributions into super- 
Hamiltonian and supermomentum are as follows: 

1 1 ■ ■ 

H_L scalar = -p^ + -99'^<^,i <^,j , (5.1) 
Hi scalar = Pfn- (5-2) 

Now, let us expand ip and its conjugate momentum p via creation and annihilation oper- 
ators: 

ip = -^2= -|- ^ -(agexp(— ifcx) + dj^exp{ikx) 

-fai" exp(zfcx) + at exp(— zfcx)), (5.3) 



P 



Po 



k>0 




—arexp{—ikx) — arexp(zfcx) 



+ai exp{ikx) + ai exp{—ikx)). 



(5.4) 
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Here V is the volume of A^-torus: 

V = {27t)^RiR2---Rn, (5.5) 



where, Ri are radii of the torus, Uj: = y 

Substituting Eqs. (|5.3|)-(p^) into Eqs. (p.lD-(pT^ and keeping only quadratic contri- 
butions one can write down the following expressions for H± scaiar{n) and Hiscaiarin) 

H± scalar n = Woifln ~ ^n) y ^ 

1 . , -+-+ ^ / V^fc^n-fc , k{n-k) 



0<fc<n \ y k n~ 



1 _ , _ I ^ ( (n + k)k ^ 

9 n+k k k n+k' I /TTi -tt: V nH 



(n + fc)A; 



fc>0 \ \/^n+kHk 



+ 2 l^K+kH + I V^^;^ + ^==== I ' (5-6) 



fc>0 



Hi scalar n ^Po (^iFi ~l~ 



0<fc<n V n-k 

k- 

+ > (a^ , ro.r + ct- , rat) — 

/ ^ \ n+fc fe n+fc fc ' 2 



fc>0 

(n + k)i 



— 'S^iara^.r + drCii ,, 

k n+k k H+k' \ ^ 

fc>0 V ^n+k 

fc>0 ^ V ^n+k 



^ n+fc K n+fc 2 V O;^ 
fc>0 V fc 

The dynamical metric does not contribute into the expansions ( p.6| ) and ( p.7| ) in quadratic 
approximation because the background value of scalar field is zero. As a matter of fact, 
such a choice of the background is the only consistent for the stationary Bianchi-J universe 
at a classical level. 

It is easy to see from Eq. (|5.6| ) that in the case when n = 0, Hamiltonian takes the 
diagonal form 

H± scalar = ^ + 3 ^ ^fc(«g«fc + ^g^fc)" (^-S) 

fc>0 

The choice of creation and annihilation operators, given by Eqs. ( ^.3|) and (|5.4| ), proves 
itself in the suggestive form of the super-Hamiltonian. 
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Using formulae (|5.6| )- ( ^.7|) one can get the expression for the contribution of the 
scalar fields into the central extension (we mean here that zero modes for scalar field and 
its conjugate momentum, as well as zero modes for other fields, are ordered according to 
Weyl rule, while other modes are Wick-ordered) 

1 Ti<i 

[L(n), ri)]c.e. scalars = '^J^[H ±{n) , Hi(—n)]c.e.scalars 



— * I A / J 1 iC\ I L I 

\n\ \ A, ^ \ ijj^ 



0<k<n 



n—k 



fc>0 



In the "last-element" limit (see sec. 3) the formula (|5.9| ) is reduced to 

1 ""-^ 1 
[L(n), L(-n)]c.e. scaiars = 2 H ^(^ ~ ^) = Y2^^^^ ~ (5.10) 

k=l 

and one can see that the contribution of scalar fields has a string-like structure. In the case 
when the model includes d massless scalar fields each of them gives the same contribution 
to the commutator [L(n), L(— n)], and thereby the right-hand side of ( p. 101 ) is multiplied 
by d. 

Consider now Wick representation for the metric and the quadratic part of gravity 
constraints. It is convenient to rewrite quadratic contribution of gravity variables into 
the constraints (|2.1| ) and (p.2|) in the following form: 



Hi = gabir'^e ,cd ,^ +9ab,e7^"'F^^\,, ,^ (5.11) 

and 

Tj iN—lry ^ab^cd , ^ r~<ab,cd,ef ^ „ i ^ T\ab,e,cdJ „ „ /'K 10^ 

H±=lp Gab,cdT^ TT + Jj^C ' ■' gab9cd,ef + ^Qab^eQcdJ- (5.12) 

Here, in the Bianchi-I flat background, we have: 

Ef,e.,^=-25;:5^)5?, (5.13) 

8cd^5^;5%^\{5%+5i5':,). (5.15) 

(jab,cd,ef _ ^ab^cd^^^ _ ^ab^{c(e^d)f _ ^{a{c ^b)d) ^ef 

_^(a{e^b)f)^cd ^ 25('^(^5^)(^5'^)/). (5.16) 

jjab,e,cd,f _ Ixabxcdr _ 1 raft r(c(e rd)/) _ 1 r(a(e red 
4 2 2 

3 1 

_|_ §(aic^b){e^d)f) _j ^{af ^b){d^c)e ^ (5-17) 



where 
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In Eqs. ( ^.13|) -( ^.17 ) we have symmetrization in couples of indices a and b, c and d and 



in Eq. ( [5 .161 ) we have also the symmetrization in indices e and /. 

Now let us expand the metric Qab and the conjugate momentum vr"^ in terms of oper- 
ators of creation and annihilation corresponding to different harmonics: 



gab{x) = y^^9abiO) + J2y^^\7^i(^ab{k)exp{-ikx) 



V ' ^ V V y2uj 

k 



k 



+a^j,(fc) exp(zA;x) + aab{k) exp{ikx) + a'^i,{k) exp{—ikx)), (5.18) 



^cd( -;\ / ^ „cd/(\\ I \ ^ -/^abjCd / ^ l^k 



= W ^i^-iy '^^^i^) + E ^g^^'^V 2/^- V V exp(-^A;f) 

VP ^ VP 

-|-a^(,(fc) exp(i/i;a;) — dab{k) exp{ikx) + a'^bi.k) exp{—ikx)) . (5.19) 



It is the choice of creation and annihilation operators which provides the diagonal form 
for the part of Hamiltonian H±[0) describing physical degrees of freedom for gravitons, 
i.e. transverse-traceless modes. 

Now one can write down the ?7,-th component of super-Hamiltonian (|4.3|) and super- 
momentum ([4.4]) in terms of expansions ( p.l8| ) and (|5.19|) : 



H^{n) = '-V2^7r''\0){-aab{n) + a+ (n)) 
^'''''^^C^''''''fgab{0)Mn) + atbin)) 



E {aab^ycd{n -k) + a+fc(fc)a+rf(n - k)) 



0<k<n 

X 



V^k^n-k ^ab,cd _|_ k)e{n fe)/ ^gfc^ed.e.f _|_ ^e^/ ^cd,ab,ef 

'^V^k^n-k "^V^k^n-k ) 

+ E (aab(A^)a+rf(n-A?) + a+fe(fc)acd(n-A?)) 

0<A:<n 

^ / V^k^n-k ^ab,cd _ (^-fc)e(^-fc)/ ^ab,cd,e/ _ fce(^ - fe)/ ^gb.ecd,/ 

+ E('^a^'(^ + k)atd{k) + a+b(^ + k)acd{k)) 



k>0 



^ / V^^H+^Qab,cd _ kekf ^ab^cd,ef _^ (^ + k)ekf ^^^^^^ ^ 

+ E('^«^'(^ + k)dcd{k) + a+b(n + k)a^^{k)) 



k>0 



V^k^n+k Qab,cd _ kekf ^ab,cd,ef _|_ (»- + fc)efc/ ^gfa^ecd,/ 
4 V^fc^n+fc V^fc^n+fc y 
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+ ^i(^abi^)(^cdin + k) + aab{k)a^^{n + k)) 

k>0 

fce(n + k)f ^ab,ecdj 



y/^k^n+k 

+ '^{aab{k)acd{n + k) + a^i^{k)a^^{n + k)) 



k>0 



4 



fc^n+fc Qab,cd _ + k)e{n + fc)/ ^ab^ed,e 



J ^ /Ce(n + k)f ^^^^^^^j 



, (5.20) 



-f^m = 5'afe(0)(a^rf(rl) - acd{n))ne'^^Ef''"^''' 



-irrcd{0){aab{n) + a+^(n)) 



'^e pab,cd,e 



+ II (a+(A?)a+^(n- A?) -a„b(A?)acd(n-fc)) 

0<A:<n 



X 



4 



T~\ab,cd,e . --c 

J^i + -r 



4 \ 



^fc _-j^ab,cd,e 



ft—k 



fee I^H-k rpab,cd,e , ('^ ^)e / "^fc rpab,cd,e 



+ XI {aab{k)a^a{n - k) - a^f^{k)acd{n - k)) 

0<k<n 



X 



^ab,cd,e ^ 



k„ [UJ. 



n—k pab,cd,e 



+ ^{a'ab{n + k)acd{k) - a^f,{n + k)a^^{k) 

k>0 

-aabin + k)a^^{k) + a+fe(n + k)acd{k)) 



X 



/ ^k rpab,cd,e + ^)e / ^fc pab.cd,. 



2 V 



fi+k 



A:>0 

-aab{k)acd{n + k) + aab{k)a^d{n + A;)) 



X 



(n + A;)e l^fi+k pab,cd,e 



n+fc pab,cd,e 



^k 



(5.21) 



Let us mention that the approximation of constraints by quadratic contribution can 
be treated as the first step of the perturbative expansion in the dimensionless parameter 



{JV-l)/2 



K 



Y{N-1)/2N- 

Indeed, due to definition of creation and annihilation operators ( ^.18 ), (|5.19 ) in which 
momentum vr*-' is proportional to 1/^p^ ^''^^ while Qij is proportional to ^''^^ the depen- 
dence on Ip disappears from super- Hamiltonian in quadratic approximation, fc-th order 
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of expansion in the metric (super-Hamiltonian contains only quadratic terms in momenta) 
is proportional to k'^"^. Thus, the terms of k-th order in creation and annihilation opera- 
tors are proportional to the same power of the expansion parameter k. So, this expansion 
in a, a^, a, a"*" can be thought about as plausible for the stationary Universe whose size is 
greater than Planck length. 

From Eqs. (|5.20|) , (|2.2| ), 1|) it follows that the contribution of graviton modes into 
the central extension of the commutator [L{n),L{—n)] looks as follows: 

[L{n) , L{—n)]c,e. gravUons = j;^ X ab,a'b'G cd,c' d' G"'^''^'^ ^ Ue'UJH 

UeUfUe' r'ab,cd,ef T?a'b',c'd',e' 

H 7^ ^ab^c'd'^cd^a'b'^-' -Tj 

ab,cdTT^a'b',c'd',e' - 



I \ " / ^ ^ab,cd 7-1' 

+ 2^ ■{-^ab,a'b'^cd,c'd'^ A 



0<k<n 
1 



, r<ab,cdrpa'b',c'd',e' , 

+ -^ab,a'b'^cd,c'd'^ '^e'^f{_k 

\^ n n r<a.b,cd,ef Ti^a'b',c'd',e'kekfkei 
+ 7:'^ab,a'b''^cd,c'd''~^ -t^i 

1^ r< r^ab,cd,ef-pa'b',c'd',e'kekf{n-k)e 

+ -^ab,c'd'^cd,a'b'^ 



, 1^ n nabM^ef ua'b',c'd',e' " k)e{n - A:)/fce' 

I 1^ n nab,cd,ef Tpa'b' ,c'd' ,e'kekfkei 
+ 7:^ab,c'd'^cd,a'b'^ 

\^ n n nab,e,cd,f j-,a'b',c'd',e' Kin - k)f{n - k)e' 

+ 7:^ab,a'b'^cd,c'd'^ -C/j 



-Mab,c'd'^cd,a'b'J^ 



, 1 ^ ^ r,ab,e,cdj T^a'b',c'd',e' K{n - k) fk^ 

+ -^ab,a'b'^cd,c'd'-l^ 



2 CJr 



, 1^ T^ab,e,cdJ Tj.a'b',c'd',e' - k)ekfke' \ 
+ -^ab,c'd'^cd,a'b'^^ -Tj > 

2 ^^fc J 

k>0 

+lGab,a'b'G,d,c'd'G'''^'''Ff''''''''''\uj^in + k),, - UJ.^^k,,) 
,1^ ^ r^gfe.cd.e/ j:.a'b',c'd',e' ( jn + k)e{n + fc)f (n + k)e' k^kfk^ 

+ -^ab,a'b'^cd,c'd'^ -C/j 







kekf{n + /c)e' 


{n + A;)e(n + fc) fk^' 




^n+k 


kekf{n + k)e' 


{n + /c)e(n + k) fke' 



< ^ n r^ab,cd,ef Tj^a'b' ,c'd' ,e' 

+ 7:^ab,a'b'^cd,c'd'^ -f i , 

2 V ^n+k 
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+ -j^^ab,c'd'^cd,a'h'J-> -C^j 



I Tiah.e.cd.f -r^a'V ,c'd' ,e' 

+ -j^^ab,a'h'^cd,c'd' J-^ 

, ^ ry r^ab,e,cd,f -r^a'V ,c'd' ,e' 

+ -^ab,c'd'^cd,a'b'^ 

f{n + k)ekfke' ke{n + k)f{n + k)^ 







^k J 


/ / 1 7 \ 7 7 

' [n + kjekfke' 


7 / 1 7 \ / 

k(,[n + k)f[n 


1 \ 

+ k)e' 




^k 




^ kf,{n + k) fke' 


{n + k)ekf{n 


+ k)e' 


V ^n+k 


^k 




^ k^{n + k) fkf.1 


(n + k)ekf{n 


+ k)e' 




^n+k 





UJ, 



(5.22) 



n+k 



It is easy to get the contractions of structure tensors G,C,D,E and F defined by 
formulae ( p.3|) , (|2.4|) , (|5.13|) - (|5.17|) . These contractions look as follows: 



r< n r^ab,cdjpa'b',c'd',e' _ iN-2){N +l) ^e' 
^ab,a'b'^cd,c'd'^ — 7^ Cj 



(5.23) 



/^ab,cd,ef rpc 
ab,a'b'^cd,c'd'^ 



cd,ef Tpa'b',c'd',e' 



ri r^ab,cd,ef rpa'V ,c'd\e' 

^ab,c'd'^cd,a'b'^ 



(5.24) 



ri n r^ab.cd.ef T^^'b' ,c'd' ,e' ^ ^ /^ab.cd.ef Tpa,'b' ,c'd' ,e' 

(Jab,a'b'(jcd,c'd'(-' -Tj — (Jab,c'd'^cd,a'b'(-' -f i 



(5.25) 



^ab,a'b'^cdfi'd'D ' ' " J^i = 7TTZ TT 0, 



A{N-1) 



4(A^- 1) 



2(iV- 1) 



(5.26) 



n ria-b,e-cdj Tpa'b',c'd',e' 3A^^ - 3A^ - 8 , f 
Gab.c'd'Gcd.a'b'-L> ' ' '■'-fcj = 7TTZ 0, d 



4(iV- 1) 



4(iV- 1) 



2(A^- 1) 



(5.27) 



Gab,a'b'Gcd,c'd'F> '^Pi = ^7777 7^ \ Oef 



1{N-1) 



Kiv-1) 



Kiv-1) 



(5.28) 



18 



^'-^.SlS^.^^!^l±^StS'". (5.29) 



8(A^-1) ' 8{N -1) 

In the "last-element limit" (see sec. 3), the expression ( p.22 ) is reduced to the much 
more simple one 

[Z/('n,),Z;( '^)]c.e. gravitons 

1 / jt— 1 oo oo 

(GGGE),,.,jin' + Y.k'-Y.k'+ E ^ 

^ \ fc=l fc=l k=n+l 

+{GGGF\,,i, (^\'pHn-k)^ 
+ {GGGE\,j, i- E f -— ^ + k{n - A;) j 

2,ttl ^+/^ A: 

+ (GGCF)i.e,. (^^' + E^') 

(71—1 OO OO 

E^' + E^'- E 
k=l k=l k=n+l 

+ (GGDF)i.e.i. (j2 kin - k)J . (5.30) 

Here the notations (GGG£')i.e.i. etc. denote the corresponding contractions taken in the 
"last-element" limit i.e. in the case when e = f = e' = N. It is easy to see that the 
coefficient at {GGGE)i,e.\. in Eq. (|5.3CI|) vanishes. The structure at {GGCE)i,e,i is very 
inconvenient for calculation, but fortunately 

{GGGE)i,^j, = 0. 

Other sums can be calculated due to simple formulae 

Efc(n-fc) = ^^^'^~^\ (5.31) 



k=l 
n-1 



n(n — l)(2n~l) 
E ^ = ^ (5-32) 



k=i 



" ,9 n(n + l)(2n + l) 

E ^ = ^ ^ -■ 5.33 

o 



k=i 
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Using Eqs. ( ^.31|) -( p.33| ) one reduces Eq. ( ^.301) to 



[-Zy(?^),Z/( ^)]c.e. gravitons 



M 3' ' 3' ' 6 

+n + i(GGCF)i.e,. + ^{GGDE\,,, - i(GGDF)i.e,.) (5.34) 

We shall be interested mainly in the coefficient at (because as it is well-known the 
coefficients at n can be adjusted to the necessary values by inclusion of intercept into Lq), 
however we shall write down the full expression for the graviton contribution into the 
quantum correction to [L{n),L[—n)\. Now we shall calculate ( [5 .341 using the following 
formulae for the contractions in the "last-element" limit: 

^PGGFU. - ^"-'f^'K (5.35) 
(GGCf ),,„,,, . J^-^)f (5.36) 

{GGDEF)u± = -i (5.37) 

(GGZ)f ),,„,,, = _M^m±^. (5.38) 



Substituting Eqs. (^35|) -(p38D into Eq. we have 



[L{n), L{-n)],,e. gravitons = + ^(^ - 2){N + 1)^ 

-J^{^ + liN-2){N + 1)^71. (5.39) 

Notice, that the appearance of the expression (A^ — 2)(A^ + l)/2 in the right-hand side of 
Eq. ( |5.39|) does not look surprising because this expression is nothing but the number of 
gravitons. What is less trivial that the coefficient at the number of gravitons differs from 
the one at the number of matter fields and even has another sign. Moreover, there is a 
constant contribution in the right-hand side of Eq. (|5.39| ), which is independent from the 
dimensionality of space. 

Concluding this section we write the contributions of some other matter fields into the 
central extension. The contribution of massless vector fields is quite analogous of that of 
scalar field. The number of local degrees of freedom of one vector field is (A^ — 1). Thus, 
if the model includes dy vector fields their contribution is 

[L{n), L{-n)l,,, sector. = ^^^^^^^n{n^ - 1). (5.40) 

The Weyl or Majorana massless spinor field has 2"^^+^^/^^^^ local degrees of freedom [] 
and the contribution of such fields into [L{n), L{—n)]c.e. looks as 

[L{n), L(-n)],.e..,„ = ^2[(^+i)/2-i]n(n^ - 1). (5.41) 



^It is interesting to recall that in the space-times of the exceptional dimensionality 2 + 8fc, with integer 
fc, one can satisfy simultaneously Weyl and Majorana conditions and the number of degrees of freedom 
of spinor is 2(^-2) for this case. For example, in the simplest 1 + 1-dimensional case (i.e. for = 1) 
the contribution of the spinor should be one-half of scalar contribution (see e.g. [p5[). 
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Let us collect the contributions of all the matter fields and gravity into the central 
extension [L{n) , L{—n)]: 

[L(n),L(-n)]c.e. = 

±n{n' - 1) + - 1) + ^2[(^+i)/2-iln(n2 - 1) 

12 ^ ' 12 ^ ' 12 ^ ' 

~Y2 + - 2)(iV + 1)) - ^ (2 + - 2)(iV + 1)) n (5.42) 

One can show by direct calculation that the commutators between area- preserving 
diffeomorphisms and between them and Virasoro-like generators do not contain quantum 
corrections. The commutators between L and L Virasoro-like generators do not include 
first quantum corrections as well. It is easy to see that the quantum correction to the 
commutator coincides with one for n)] given by the right-hand 

side of Eq. 

Thus we have seen that the algebra of constraints is disclosed already by constant 
contributions of the first quantum corrections. It is worth mentioning that the fact of 
a disclosure does not mean by itself the quantum inconsistency of the theory. A proper 
test for consistency is to examine the nilpotency condition for the quantum BRST charge 
related to the disclosed constraint algebra. This examination is performed in the next 
section. 



6 BRST operator and critical relation 
between parameters of the theory 

In this section we study the operatorial canonical BFV quantization [|, ^ of the con- 
strained dynamics developed in the previous section for the toroidal Bianchi-I cosmology. 
The general prescription of the BFV-method implies making a definite choice of the op- 
erator symbol both for the original phase space variables and ghosts. The main insight 
at the choice of the ghost ordering is from the due regard to the structure of the related 
constraints. The Virasoro-like subset of the constraints has the polarization, splitting 
them into the positive- and negative-frequency ones, being conjugated to each other. 
Thereby, it is natural to subject the related ghost variables to the Wick ordering rule. 
Meanwhile the ghosts attached to the area-preserving diffeomorphisms are quantized by 
the Weyl rule. There are several reasons for such a choice of ordering for these ghosts. On 
the one hand, this ordering, as it will be shown later, is compatible with the mentioned 
in previous section fact, that there are no quantum corrections to the commutators of 
the area-preserving diffeomorphisms among themselves and with the Virasoro-like con- 
straints. On the other hand, this choice follows to some extent to a certain analogy with 
p-brane theory where the area preserving constraints describe some residual sym- 
metry after fixing of the light-cone gauge. Now we are in a position to write down the 
quantum BRST charge for the model. It looks like 

= ^(c+(")L(fi) + c(")L(-fr) + c(")Z(n) + c+(")L(-n)) 

n>0 
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n,v 

n>0,m>0 n>0,m>0 

+ 2 ^ Mn+m,^')^L(n)L(m) + o ^ M"+H^L(n)Z(m) 

n>0,m>0,i; n<0,m<0 

^2 -^(n+m)^L(n)Z,(m) ^ 2 ^ -^("+'^.'!')^L(ri)L(?n) 

n<0,m<0 n<0,?n<0,i; 

^ 2^ ^ ^{n-m)C I-' L{n)L{-m) ^ ^ ^{n-m)'- L(n)L{-m) 

n>0,n>rh>0 n>0,n>rh>0 

+7 V r+('^V(*)p+ TjH{fi-m,v) .,?=("^)7-/^("-"^) 

^{n-m,v)'-' L{n)L{-m) ^ ^ ^{n-m)C 'J L(n)L(-m.) 

fL>Q .n>fh>{) ,v n<0,n<m<0 

I V r+i'^)p+ r^m)TjL{n-_m) ■ +{fi)^{m)p+ jjH(n-m,v) 

^ 2^ ^ ^(n-m)^ ^L(n)L{-rn) ^ * C C -f^(n-m,'iT)^L(n)L(-m) 

n<0,n<rn<0 n<0,n<m<0,S' 

n>0 n>0 

n<0 n<0 

^ ^in+m)^^ ^L{n)L{m) ^ 2^ ^ ^{n+m)^ ^ L{n)L{m) 

n>0,m>0 n>0,m>0 

I,- V r+^'^^P/-, - -xr('^)r/'^^"+'^''^^ + V r+('^)r-+(") P.- jMn-m) 
-rt ^ C -f^(n+m,'y)t' ^L(n)L(m) ^ Z^ -^("-H^L(n)L(-m) 

n>0,m>0,'!/ n>0,n>m>0 

V 7^+(™)r+(")p+ TjUf^-m) ■ sr^ +(rn:)+(n)p ^ jjH{n--m,v) 

-r 2^ C L. ^{n-m)'-' L{n)L{-m) ^ * Z^ M"-"^,■f)^L(n)L(-m) 
n>0,n>m>0 n>0,n>rn>0,v 

/ (n—m) L{n)L{—m) Z-^ ""-J L(n)L{—m) 

rh>0,rh>n>0 rn>0,m>n>0 

+7 V rH-(m)^+{n) 6_ ^ ^/-/'^('^-™'^) I V r+^^^r^"^'") P.- ^,7'r^("+™) 

_^ Z^^ ^ -^("-m,'y)t^2.(n)L(-m) ^ ^ 2-^ ^ ('^+'^) L{n)H(rn,v) 

m>0,m>fi>0,v n>0,m>—n 

-^ V P+ ^+(")Mm.i?)7-r^("+"^) _L,- V r-+("),^('^''^)p.- - TrH{n+m,w) 

n>0,m.>— n n>0,m>— n 



n+m) 
)^(m,i;) 



+ ^ Z^ / ^ Mn+H^L(n)^(m,i;) + ^ Z^ /(n+m)^ C U 
n<0,m>n n<0,rh>n 

n<0,m>n n>0 

n<0 n>0,7n<— n 

yn-tm) L{n)H{rn,v) ^ yn^ra) L(^f{^H{m,v) 

n>0,m<— n n>0,m<—n 

-^ V P+ ^('^i>'^)7-7"^("+"^) -U V r+in)p;{m,v)p ^TTL{n+m) 

^ ^ 2^ f{n+ih)^ ^ ^L{n)H{m,v)^ ^ 4^ ^ M"+"»)^L(n)^(m,i;) 
n<0,m>— n n<0,m>—n 
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n<0,m>— n 

+h.c. 

\^ ( pi{rn,w) ~(m,w) p ^ ^, _|_ P,^ ^ . 

+e<•->A«+«,.^c'--^)4'?;■Sfi,,,. (6.1) 

In Eq. ( |6.1| ) we should make summation over the set of (A^ — 1) unit vectors v enumerating 
area-preserving diffeomorphisms constraints H^P^ which are orthogonal to vector n. 

The ghosts c'^^'^\ c^^'^\ P^t^ and P^-^ are creation operators, c^"-*, c^"\ P^it) and P(^) are 

annihilation operators. All these operators are Wick-ordered, while c^'^''"\ P(^^^), c^^\ P(o), c^^^ 
and P(o) are ordered according to Weyl rule. The latter are subjected the following non- 
vanishing anticommutation relations: 



C , -r(m)J — O(^) 



P+ 1 



^(") P+ 1 



(m)' 



(6.2) 



c°,Po]=^, 

:c°,Po] = ^. 



Making use of Wick theorem, one may establish that due to the terms with two con- 
tractions of ghost operators, the requirement of vanishing of the squared BRST operator 

= should give rise to the following nonvanishing corrections to the quantum involu- 
tion relations of the constraint operators: 



[L{n),L{-n)]c.e. ghosts 
^ — ' Lin 



2 L{n) L{Q)^ L{n) L{~n) 2 H^^) L{-n)'^ L{0) L{-n) 



0<k<n 



L{n-k) TjL{-k) 
L{n) L(-k) L{n~k) L{-n) 



L(n) L(-k 

_i_Sr^(TrL{n+k) jjL{k) jj-ny-Kj jj 

L{n) L{k) L{n+k) L{-n)^ ^ ^ ' ^ 



A:>0 



L(—k) jjL{—n—k) ^ 

L{n) L{-n-k) L{~k) L{~n)' 



(6.3) 



[I^y I L l^yii)\c.e. ghosts — 2^ U-n) L{0)^ L{-n) L{n) ^ 2 ^(-") -^^W 

^ L{-n)L{~k) L{-n-k) L{n) 
0<k>~n 

,Sr^(TM-n+k) TjL{k) ^ iTT^i^) _ TjL{"n~k) n ((\ A\ 

~^ Z^'^^L(-n) L(k) L{-n+k) L{n) L(-n) L{n-k) L(-k) L{-n)>' ^ ' 

k<0 



23 



Straightforward but tedious consideration shows that the nonvanishing quantum cor- 
rections arise only in commutators [L{n), L{—n)] and [L{—n), L{n)]. Thus, the rest of the 
commutators do not get c-number quantum corrections. Moreover, all the contributions 
to Eqs. ( |6.3| ), ( |6.4| ) of structure constants involving area- preserving diffeomorphisms are 
mutually canceled. 

Left-hand sides of the relations (|6.3|) and (6^) represent the quantum corrections to 



the corresponding commutators calculated in the preceding section, while the right-hand 
sides consist of the quantum ghost contributions required by the nilpotency condition. 
^ Left-hand sides are determined by a specific realization of the constraint algebra and, 
in particular, depend upon the specific spectrum of the fields involved into the original 
Lagrangian, whereas the right-hand sides do not depend on the specific realization of the 
constraints, being completely fixed only by the structure constants of the algebra. So, 
these equalities, following the nilpotency of the quantum BRST charge, may constitute 
in principle the nontrivial restrictions on the theory parameters. 

It was shown in the preceding section non-trivial corrections to quantum involution 
relations (see Eq. ( ^.42|) depending on matter content of the theory appear only in 
commutators [L(n), L(— n)] (and, symmetrically, in [L{—n),L{n)]). Here, we have an 
analogous structure of quantum corrections to involution relations being originated from 
quantum commutators between the terms in Q which are of third order in ghost variables. 
In this approximation, it is sufficient to keep only the constant parts of the structure 
coefficients because the multipole harmonics may contribute into the central extension 
only if the terms with three or more contractions are accounted for. These terms represent 
the quantum corrections of a higher order. 

Now we should calculate the right-hand side of Eq. ( p. 3D explicitly (calculation of 
( |6.4| ) does not give us an additional information). 

Substituting into Eq. (|6.3|) explicit expressions for structure constant we have: 



[L(n), L(-n)],.e. ghosts = 2n^ + ^ ..^.r-,}^ — r^ 

X |2|n — k\ff'k'^ — \n — k\n^k^ — \n — k\n^k'^[nk) 

+2\n - k\ff{nkf + 2\n - k\k^{nk) 
-3|n - k\P{nkf -\n- k\{nkf + An'^\k\{nk) 
+n^\k\^ - n^\k\^{nk) - Qn^\k\{nkf 
-2\k\^{nk) + h\k\^{nk:f - \k\{nkf] 

^^2\n\k\n + kY 

X |2|n + k\n^}? -\n + k\ffk'^ + \n + k\ffl?{nk) 

+2\n + k\ff{nkf - 2\n + k\k'^{nk) 

-3\n + k\P{nkf + \n + k\{nkf - An^\k\{nk) 

+ff\k\^ + fi^\k\^{nk) - Qn^\k\{nkY 



2 



2r 



The general fact of appearance the corrections which emerge from the quantum contribution of ghost 
commutators has been first studied in the paper j?] . 
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+2\kWnk) + 5\k\%nky + \k\{nkf} 



(6.5) 



Going in Eq. (|6.5|) to the "last-element limit" we come to the following expression 

[L{n),L{-n)]c.e. ghosts = + ^(2ra - k){n + k) = —n^ - -n. (6.6) 

k=i " " 

It is easy to see that the result ( |6.6|) coincides with the well-known expression for the 
Virasoro algebra. 

Now, equating the coefficient at with an analogous coefficient in Eq. ( |5.42|) we 
come to the following critical relation between the dimensionality of space and the matter 
content of the theory: 

d + dv{N - 1) + dp2[(^-i)/2-i] _ 30 + ^^jv -2){N + 1). (6.7) 

Thus, we may claim that this short relation eventually represents a necessary condition 
for a nilpotency of the BRST charge ( |6.1|) in the first quantum approximation. As is seen 
it imposes some correlations between number of scalars d, vectors dy, spinors dp and 
dimensionality of the space. 



7 Conclusions 

In the Conclusions let us summarize briefly the key points of suggested construction and 
make some comments about the results. 

We have suggested a treatment of closed cosmological models by the harmonic expan- 
sion of gravity and matter canonical variables and the constraints. As basis harmonics 
we have chosen eigenfunctions of Laplacian on the maximally symmetric manifold of the 
topology under consideration. In doing so the structure constants of the involution re- 
lations of constraints algebra could be expressed via Clebsch- Gordan coefficients of the 
corresponding symmetry group. This provides the general background for the multipole 
expansion of constrained Hamiltonian dynamics of closed cosmological models, and their 
subsequent quantization. In this paper, however, we restricted the further analysis to the 
case of - dimensional generalization of the stationary closed Bianchi - / cosmological 
model (A^ - torus). Except the technical simplicity, this choice of the particular model 
was motivated by the fact that it has the stationary classical solutions. 

Besides the multipole expansion, the construction implies splitting of the constraint 
set into the area-preserving diffeomorphism generators, forming a closed subalgebra, and 
the Virasoro-like generators. 

Performing the quantization we have constructed the BRST operator, defining the 
Wick ordering for the gravity, matter and Virasoro-like ghost excitations and the Weyl 
one for zero-modes and for the ghosts related to area-preserving diffeomorphisms. Our 
further objective was an examination of the nilpotency condition for the BRST charge. 
Actually, in the nilpotency equation = 0, we have managed to calculate the constant 
part of the first quantum correction to the coefficient quadratic in ghosts . In doing so 
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we applied the expansion of constraints and structure functions in powers of the dimen- 

;(JV-l)/2 

sionless parameter k = y(^,i)/2jv • Finally, we derived from the nilpotency condition the 
relation between the dimensionality of the space and the spectrum of matter fields 
(|6.7|). Although, the restrictions imposed by Eq. (|6.7| ) do not predict the unique spec- 
trum and spatial dimensionality, one can observe some curious consequences, for example: 
stationary Bianchi-/ closed cosmology is inconsistent without matter at any positive in- 
teger N. Notice, also, that in = 3 case, the admissible matter spectra do not contain 
sufficient number of degrees of freedom to be compatible with the Standard Model or 
its generalizations. It is not too disappointing, because the stationary closed Bianchi-/ 
model itself could hardly be expected to describe the observable cosmology. On the other 
hand, the very existence of such a relation allows to hope that analogous conditions, being 
deduced for more realistic models (with non-stationarity, other topology, massive particles 
and, perhaps, with an extended symmetry), can display a better correspondence with the 
reahstic particle physics. 

One more (a little bit speculative) way of fine-tuning the number of matter degrees 
of freedom to the dimensionality of space could be found in a generalized Kaluza-Klein 
ideology. In this way, Eq. ( |6.7D can be thought of as valid for the fundamental multi- 
dimensional theory. That provides much more opportunities for the matter spectrum in 
low-dimensional effective theory by means of a proper compactification. 

While an opportunity of interplay between matter content of the Universe and its 
geometrical characteristics is of certain interest, we would like to think that the designed 
scheme of Hamiltonian treatment and quantization of closed cosmological models may 
open some other promising prospects for studying various aspects of quantum cosmology. 
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